Theory of superconductor with k close to 1/v^ 
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As was firstly shown by E. Bogomolny, the critical Ginzburg-Landay (GL) parameter k — l/\/2 
at which a superconductor changes its behavior from type-I to type-II, is the special highly degen- 
erate point where Abrikosov vortices do not interact and therefore all vortex states have the same 
energy. Developing a secular perturbation theory we studied how this degeneracy is lifted when k 
is slightly different from l/-\/2 or when the GL theory is extended to the higher in T — Tc terms. 
We constructed a simple secular functional, that depends only on few experimentally measurable 
phenomenological parameters and therefore is quite efficient to study the vortex state of supercon- 
ductor in this transitional region of n. Basing on this, we calculated such vortex state properties 
as: critical fields, energy of the normal-superconductor interface, energy of the vortex lattice, vortex 
interaction energy etc. and compared them with previous results that were based on bulky solutions 
of GL equations. 



I. INTRODUCTION 

Although the Ginzburg-Landau (GL) theory covers 
all the variety of superconductors, both of the type I 
with GL parameter n < l/\/2 and of the type II with 
K > l/\/2, the most theoretical studies of the vortex 
state deal with the case of k ^ l/\/2 since at k ^ l/\/2 
the GL equations are simplified substantially and also 
the grand majority of type-II superconducting materials, 
including High-Tc superconductors, correspond to this 
limit. 

Studies of superconductors in the transitional region of 
small K ~ l/\/2 done mostly in the 70-s were based on the 
solution of the full system of GL equations and require 
bulky calculations. Meanwhile, E. Bogomolny proposed 
in 1976 an elegant jWay to operate with similar problem 
of the string theorya and showed that at the special point 
K = l/\/2 the order of the equations can be reduced and 
all the vortex states with arbitrary located vortices have 
the same energy when the applied field is equal to the 
critical field He- 

Historically the Bogomolny approach was done for the 
high-energy physics and the superconducting commuiiity 
was unaware of it even when L. Jacobs and C. Rebbin re- 
formulated the Bogomolny equations in terms of super- 
conductivity and demonstrated that they can be writ- 
ten in a form of nonlinear electrostatic equations of the 
Boltzman plasma that we are calling as Bogomolny Ja- 
cobs and Rebbi (BJR) equation. Only very recently the 
Bogomolny Biethod was used to study vortices in meso- 
scopic disksB and to calculate the structure of multi- 
quanta vorticesu in superconductors with k — l/\/2- 

In the present paper we derive a regular way to treat 
the vortex state of superconductors with k close to 
l/-\/2 basing on the BJR approach. Our basic idea is 
to consider the highly degenerate Bogomolny state at 
K = I/V2, H = He as zero approximation and then 



to account for deviation from this point via the secular 
perturbation method that lifts the degeneracy and selects 
the most stable vortex configuration. In Sec. Ill we give 
a brief overview of the Bogomolny method and then, in 
Sec. [V, present the solutions of BJR equation for the 
different vortex lattices. In Sec. we develop the secular 
perturbation approach. The following perturbations can 
lift the degeneracy. 

(i) Deviation of k from l/\/2 that is accounted by a 
small parameter. 

7==k2- - « V2k- 1. (1) 

(ii) Deviation of the applied field from He- 
(iii) Next in 

t = T/Tc - 1 (2) 

corrections to the GL functional. 

(iy) Thermal fluctuation effects. 

(y) Finite size and demagnetization effects. 

We consider only the first three contributions and show 
that they can be incorporated in a very simple secular 
functional that acts on the Bogomolny degenerate solu- 
tions and looks like a six-order polynomial for the am- 
plitude of the order parameter. This functional depends 
only on few phenomenological parameters that can be 
found from experiment and that completely determine 
the behavior of superconductors with k ^ l/\/2 in a 
magnetic field. 

Certain properties of superconductors with k ^ 1/V2 
were calculated either from the GL theory extended to 
low temperatures or from the microscopic Gorkov equa- 
tions. These calculations, overviewed in Sec. O, were 
however dealing either with cumbersome analytical ex- 
pansions or with numerical computations that both are 
difficult to catch on. It is therefore of interest to recal- 
culate these properties in a systematic perturbation way 
and compare them with the older results. 



In Sec. Vl we calculate the following parameters of a 
superconductor with k ~ l/\/2: 

a) Critical fields Hd, i?c2 and He- 

b) Energy of the NS interface. 

c) Energy of the regular vortex lattice as function of 
the applied field. 

d) Energy of the N-quanta vortex. 

e) Vortex interaction that can have unconventional at- 
tractive character. 

Basing on these calculations we discuss the possible 
scenarios of the normal-superconducting (N-S) transition 
in a mag netic field for a superconductor with k ^ l/\/2 
(Sec. [VIl| ) that occurs either directly (like in a type-I su- 
perconductor) or via formation of the intermediate vor- 
tex (V) state (like in a type-II superconductor). The 
actual scenario depends on the relative strength and sign 
of coefficients in the perturbation functional that can be 
extracted from experiment. We calculate the location 
of the triple point L on the H — T plane where the di- 
rect N-S transition splits into N-V and V-S transitions 
and a superconductor changes its behavior from type-I 
to type-II. The important feature is that, the N-V and 
V-S transitions close to point L can be both continuous 
and discontinuous unlike the traditional type-II super- 
conductor with K ^ l/"\/2. We calculated the location 
of tricritical points T2 and Ti where the N-V and V- 
S transitions change their character from continuous to 
discontinuous. 



II. PREVIOUS STUDY 



A. Theory 

Already in his pioneering workS Abrikosov noted that 
solution of the GL equations at k ~ l/v2 is a sepa- 
rate and quite complicate problem. Since that various 
related theoretical investigations that are partially re- 
viewed inETS were done. The first series of investigations 
dealt with an expansion of the BCS free energy close 
to Hc2 over a small parameter H — Hc2- The magnetic 
and thermodynamic properties of-superconductoij-close 
to Hc2 were calculated for dirtyD'EEl and intrinsicO su- 
perconductors. The most complete calculations of this 



type are given inlia. The possibility to have a discontin- 
uous N-V transition iuj-a superconductor with k ^ l/\/2 
was firstly indicated irJa. 

On the basis of BCS theory Tewordt and Neu- 
mann calculated the low-temperatuife,correcti(ms to GL 
functionalljilj and found the upperE£l and IowEj critical 
fields with an accuracy t^ at arbitrary k. j— . 

Basing on this extended GL functional, A. E. JacobstZl 
considered a superconductor with k ~ 1/V2 and calcu- 
lated the NS interface energy, the energy of single and 
double quantized vortex. He obtained that at certain 
conditions the vortices in a type-II superconductor at- 
tract each other and predicted the discontinuity of the 



V-S and N-V transiiipns. The analogous result was also 
obtained by HubertES . 

Grofimann and WisseO calculated the free energy of 
a superconductor with k ^ l/v2 close to Hc2 using the 
extended (although not complete) functional of Tewordt 
and Neumann. They founded a discontinuity of V-S tran- 
sition in a limit of the dense vortex lattice. All the above 
conclusions were reproduced by BrandtnJ who developed 
a variational numerical method to solve the Gorkov's 
equation for vortex lattices for all possible values of H, 
T and k. 

Recently OvchinnikovEJ carefully derived the coeffi- 
cients of the extended GL functional from a microscopic 
theory for different types of the electron scattering. He 
considered an expansion of the free energy near Hc2 up 
to the order of {H — Hc2Y and specified the case when at 
K ^ l/\/2 the N-V transition has a discontinuous char- 
acter. 

In the present article we reproduce the above results in 
a more simple way, basing on the Bogomolny treatment 
of superconductors with n = l/\/2. 



B. Experiment 



The superconducting metals Ta, Nb, In and Pb with 
K close to l/\/2 were intensively studied in the 60-s and 
70-s. The variation of k was achieved either by dissolv- 
ing of foreign atoms of N, Tl, Bi or by preparation of 
samples with, different defect concentration. We 
to magneti(£3 calorimetricc3 and neutron diffractioncJ 
experiments in pure Nb (k ~ 0.85 — 0.96),|-Jto mag- 
netic measurements in TaN (k ~ 0.35 — 1.53)123, in Nb 
(k ~ 0.78- 1.03)Ea and in InBi (k ~ 0.76- 1.46)kzl and to. 
direct observation of vortices ip|EJpTl (k ~ 0.43 — 1.04)l3 
and in Pbin (k '--^ 0.76 — 1.46)ol2J by decoration. Jiefer- 
ences to other related experiments can be fond iiol3. 

The fact that the V-S transition can be of the first 
order at k ~ l/-\/2 was discovered already |-iji—jliie 
early magnetic and thermodynamic experimentfS'OtZl. 
The detailed magnetic study of a superconductor that 
changes its behavior from type I to type II was done 
for tantaljam samples with some amount of dissolved 
nitrogenc3. A discontinuity of the vortex lattice pa- 
rameter at the V-S transition was observed in neutron- 
scattering experimentgHilS 

The convincing confirmation of discontinuity of the V- 
S transition in superconductors with k ^ l/v2 was done 
by a direct obser|5iation the vortex domains inside the 
Meissner phasecfl Lj. Such coexistence of different phases 
is known to be a signature of the first order transition be- 
tween them. This Jatermediate-mixed domain structure 
was interpreted inc3 in terms of a long-range vortex at- 
traction. 

The discontinuity of the V-S transition provided by 
an attractive interaction between vortices is therefore a 
well established fact. Meanwhile the ground state of 



the vortex lattice and the configuration of domains of 
the mixed-intermediate phase jajre-still unclear. Although 
the decoration experimenta3Cl^c2rE3 demonstrate the very 
peculiar magnetic textures, including vortex segregation 
and clustering into lamellar and drop-like domains, no 
systematic study of this question that take into account 
the demagnetization and finite-size effects was done. We 
believe that our calculations of the vortex energy in the 
bulk superconductor with k ~ l/-\/2 can be extended to 
simulation of magnetic textures in the realistic finite-size 
samples. 



III. GL FUNCTIONAL AT k^ = 1/2 AND BJR 
EQUATION 

In this Section we describe the Bogomolny proceduretl 
that allows to simplify the GL functional and to reduce 
the or|der of the GL equations at k = l/v^. Jacobs and 
Rebbitl formulated the Bogomolny equations in a simple 
form of the nonlinear Poisson equation that we shall call 
as BJR equation. We discuss the properties of the vortex 
solutions of the BJR equation and their interpretation in 
terms of electrons in Boltzman plasma given incl. 

We start from the conventional GL functional: 



j^^am' + ^m* + K\B^\' + — -^^' 



8tt 



Att 



(3) 



where 



2e 
a = ait, D = V — i — A, B — rotA. 
ch 



Refer first to the characteristic parameters of supercon- 
ductor. In the uniform state the superconducting order 
parameter takes the equilibrium value: 



1/2 



Ratio of the penetration depth and coherence length: 



S= - 



32TrKe^ a 
gives the GL parameter: 
S 

K = — = 



a 



chg^/^ 



e (327r)i/2 \e\K ' 



(5) 



(6) 



Thermodynamic critical field and the upper critical field 
are written as: 



Hc = 



ch 



— a, Hc2^- ^. I „ < 

9 J 2|e|L> 



= V2kHc 



(7) 



Note that the commonly used expression for the low crit- 
ical field: 



Mci — -: — ;r77 In /x — Jlc 



AttS^ 



V2« 



(8) 



is valid for k ^ l/-\/2 and is not applicable in our case. 
Corresponding expression for Hd at k ~ 1/a/2 will be 
obtained in Sec. V|. Introduce now the dimensionless 
variables that are slightly different from the commonly 
used in the GL theory: 



TD 

b = V2n — , 



R 



H, 



kA 



hr, = V Zk — , a = - — 



H, 



5Hr 



/ = 



T 



2 , 2 



The GL functional O) in this variables takes the form 



(9) 



/ = n'M'' - 1)' + |( V - za)^|2 + (y - bh^). 



(10) 



It is convenient to use the complex variables: 



1 



(11) 



a,9=-(V,TJVj,), 



a,a= ^(oa; Tiay) 



in which the GL functional is written as: 
/= 2\{d - ia)ij\^ + 2\(d~ia)tp\ 



wm 



I) 



9 



(12) 
(13) 



To catch the special properties of the GL functional at 
n = l/\/2 one can integrate the first term in ( p^ ) by 
parts. Using that: 



j fdgdS = - j 



gdfdS+'-jfgdC 



^ fdgdS 

and that: 

b = —2i{da — da) 
one gets the substitution: 



2\{d-ia)tP\^ ^ 2 Ud- 1^)1^1^ + \ij\^b. 



(14) 



(15) 



(16) 



We neglect the contribution of the surface currents that 
are important for finite-size effects considered irjj. Fi- 
nally, one comes to the alternative expression for /: 



+ 7(IV^l'-l)' +il-ho)b, 



(17) 



where 7 = k^ — 1/2. When 7 = and /iq = 1 (i-e- 
H = He) the functional (17) reduces to the sum of two 
square terms. The absolute minimum is achieved when 
both these terms are equal to zero i.e. when the following 
equations are satisfied: 



{d ~ ia)ijj — 



and 



(18) 



(19) 



1 — ^ip = —2i{da — da) — b. 
Substitution of a from (jl^) : 

a^-id\nip (20) 

to (n9) gives the BJR equation: 



iv^lnlV^r 



\i:\''-l + 2nY,N^S{r-ri). 



(21) 



Firstly introduced inQ the (5-function terms correspond 
to the iVi-quanta vortices located at r = r^ where ■0 gets 
the phase wind 2TrNi and In |V^| has a In r^ singularity. As 
follows from Eq. ( [l9|) the distribution of magnetic field 
inside a sample is uniquely related with the amplitude of 
the order parameter: 

6(r) = l-|^(r)^ (22) 

The BJR equation can be alternatively written as: 

V^(^ = e2'^-l + 27r^7Vi5(r-r,). (23) 

where (p — In \ip\. This form has a simple interpretationQ. 
It describes the screening of electrons in a classical Boltz- 
man plasma that consists of the positive ionic background 
with potential l/47r. 

BJR equation ( pi|) defines the amplitude of the super- 
conducting order parameter |'(/'(r)| and the magnetic field 
6(r) at 7 = and at ft-o = 1 as a function of position 
of the vortices ri...r„. All these vortex solutions corre- 
spond to the absolute minimum of the functional ( |l7| ) 
and, therefore, 7 = and /iq = 1 is a special highly de- 
generate point where all the vortex states have the same 
energy. 

This infinite degeneracy over ri...r„ is lifted if one goes 
either beyond 7 = 0, /iq = 1 or beyond the GL approxi- 
mation. When 7 = and ho > 1 the absolute minimum 
of ( [l7| ) is the normal state with b — h^ and ji/'l — 0. When 
7 = and /iq < 1 the absolute minimum of (|l^ is the 
uniform superconducting state with 6 = and \^\ = 1. 
To find the vortex states when 7 f^ one should account 
for the term ^{\'4>\ — 1)^ as the secular perturbation that 
lifts the degeneracy. This will be done in Sec. ^ together 
with account of the low temperature corrections to the 
GL functional. 



IV. VORTEX STATE AT k^ = 1/2 

A. General 

In this Section we discuss a particular class of solu- 
tions of BJR equation ( pl| ) when N— quanta vortices are 
packed into the regular lattice with basis vectors ai,a2. 
We will need these solutions in Sec. VI as zero approxi- 



mation of the perturbation theory to find the most stable 
vortex configuration beyond the Bogomolny point. The 
unit cell area S = aia2 sin a {a = ai'^a2) carries the fiux 
27rA^ and therefore is related with the average induction 
as: 



S = 2-KN/b. 



(24) 



The value of b and S varies from 6 = 0, S* = cx) (al- 
most non-overlapping vortices) to 6 = 1, ^ = 2tiN (dense 
vortex lattice). We consider both the limits analytically. 
We use a numerical procedure to treat the case of an 
arbitrary lattice. 



B. One vortex 

The axially-symmetric distribution of the order param- 
eter gAr(r) = |'0(^)| inside the iV-quanta vortex is calcu- 
lated from the radial version of BJR equation (21): 



II 


9'?i 


9'n , 3 


9n - 




+ 9n 




9n 


r 



9N 



(25) 



with the boundary conditions: gAf(O) = 0, gN{oo) = 1 
and asymptotes: 



9N 

9n 



Bnt 



\-A 



= -V2r 



N- 



' ^ 0, (26) 

-^00. (27) 

1 and got Bi sa 0.9 



We solved (Eq) numerically for N 
and Ai w 1.6. 

To find the vortex solution at iV > 1 it is more con- 
venient to use the new function W7v(r) = gjj (r) that 
satisfies the equation: 



'^N 



12 



VN 



'-"N 



1 

TV 



,2N+1 



1 

N 



VN 



(28) 



and has a linear behavior ^ i?^ r at ?' — > 0. The an 
lytical expression for gNW) at TV ^ 1 was obtained i: 
In dimensionless units (M it is written as: 



9N = \ — 

JN 



'i(-'-^) 



r < TAT. 



The size of the vortex core 

rN ^ 



(29) 



(30) 



is estimated from that, the almost uniform magnetic field 
the flux 27rA^ 



1 , distributed inside the vortex area nrj^ results to 



C. Separated vortices and diluted lattice 

The magnetic flux of slightly overlapping vortices can 
be written as the superposition of fluxes of separate vor- 
tices: 

6(r) = ^6^(r-rO = ^(l-5^W), 

i (31) 

where gN{i) = 5iv(|r — r^j ) is the solution of (p5|). Then, 
the amplitude of the order parameter is written as: 



|^(r, v,...v^)\^ = 1 - 6(r) = 1 + ^ {glii) ~ l) 



(32) 



D. Vortex bunch 



Consider now the group of A^ vortices located close to 
the origin such that |ri| <C 1. This vortex bunch can 
be viewed as the TV— quanta vortex (7Ar(r) with the split 
core. By direct substitution one proves that the corre- 
sponding solution of equation (p3) within the accuracy 
0(max Ir^l ) is given by: 



(/7(r,ri...rAr) = — ^^jv(r-r,) (33) 

i 



N 



W^lnlV-ol 



1 



2^Y.^{i 



(37) 



that close to Hc2 coincides with BJR equation in a limit 
< 1-01 >^ 0, 6 — > 1. To find the normalization coef- 
ficient A(b) ^ (1 — 6)^/^ one should treat the nonlinear 
part of BJR equation as a perturbation. 

The A^-quanta lattice solution with unit cell area 2-kN 
can be written in the analogous way: 



V'o(r) = AN{h)e'\zdhT"/2TTN,T)e 



=-br/2 



(38) 



F. Arbitrary lattice (numerical solution) 

We performed the numerical integration of BJR equa- 
tion for square and triangular vortex lattices with N — 
1 , 2 in a whole interval of < 6 < 1 presented it in a more 
suitable form. First we pick the zeroes of the order pa- 
rameter via the special multiplier ip^ (r) that was taken as 
( p5| ) for iV = 1 or as (|3^) for an arbitrary N and present 
the order parameter in the form: 'ip{v) = |V'o(r)| • |'!/''(r)|- 
The new equation for function |V''(r)| has no singular 
i5-function terms and is written as: 

iv^lnlV-f = |0o|'|V'f -(1-6). (39) 



Taking ip' — In \ip'\ we present (39) in a form 



where ^N{r) = ln|gjv(?')|, <^(r, ri...r7v) == 
In |-0(r, ri...r7v)|; the origin being taken in the center 
of the vortex gravity such that Eir^ = 0. 

With the same accuracy the order parameter is written 
as: 



(r, ri...rAr)| = g^ir) + :rT7ffw(0 V(riV)2 ln5w(r). 

^^ r (34) 



E. Dense lattice 

The order parameter of the dense 1-quanta vortex lat- 
tice is presented by the Abrikosov solution close to Hc2'- 



V'o(r) = A{b)e{z^JhT"/2TT,T)e-''y'l^, 



(35) 



where: t = t' + it" — 026*" /oi is the geometrical pa- 
rameter of the lattice cell (for square lattice r = i; for 
triangular lattice: r = e^'"/'^) and 9 is the Jacobi theta- 
function: 

00 
0(z,t) = 2^(-l)"e*""("+i/2)%ij^[^(2n + l)z]. 

n=0 (36) 

Function 0o(r) satisfies the linear equation: 



V2^' = |0o|'e''^'-(l-fo). 



(40) 



This nonlinear Poisson-like equation can have the peri- 
odic solution only if the electroneutrality condition is sat- 
isfied: 



1 



=<|^o|'e'^' 



> . 



(41) 



Taking into account ( p4| ) and performing the rescaling 
r ^ r • (27rAf/5)i/2 we map Eq. (E^) onto 



VV = A ( 



lV'o(r)| 



^eV 



< \My)\ e^'^' > 



1) 



(42) 



that is defined for the parallelogram of the fixed area 
S = 1 with periodic boundary conditions. The parame- 
ter of the equation 



A = 27rA^ ( = - 1 
b 



(43) 



varies from at Hc2 to 00 at Hd- 

The problem was solved for the square and triangu- 
lar vortex lattices with A^ = 1, 2 using the Matlab PDF 
Toolbox by the finite element method with the adaptive 
mesh refinement and with the rapidly converging Gauss- 
Newton iterations that were used to account the non- 
linear right-hand-aide of (^4). For details of the numerical 
method see RefO. The obtained solutions were verified 



by the direct substitution them back to (42) 



G. Normal-superconducting interface 

The profile of the NS interface is usuaUy considered in 
two limit cases: when k — > oo or when k — > 0. It ap- 
pears, however, that the NS profile can be found exactly 
at K = l/v^ by integration of the equation (^3|) that in 
ID case looks like: 



" = e^^ - 1. 



The first integral of (H): 
alternatively can be written as: 



(44) 
(45) 
(46) 



The integration constant in (Eq), (Eoh was chosen to sat- 
isfy the NS interface boundary conditions: 

IV-I = 0, d\ip\/dx^O whcnx^-oo, (47) 

{ipl — 1, dlipl/dx — O when a; ^ cxD. 

Further integration of (Bq) gives the implicit form of 

[■(/'(a:^)! at the NS interface: 



IV-I 



dy 



Y/t/"* - 2/2(1 + In j/2) 



(48) 



V. PERTURBATION THEORY 

To find the most stable configuration of vortices be- 
yond the infinitely degenerate point 7 = 0, ft-o = 1 of 
Bogomolny functional dl^) we construct the secular per- 
turbation functional that acts on the (zero order) de- 
generate solutions \^{r,ri...rn)\ of ( pl[ ) and selects the 
vortex configuration ri...r„ having the lowest energy. 

The perturbation for 7 and h were given already by 
two last terms in wTvj. To find the perturbation for t 
one should ext end , ±he GL functional to low tepipera- 
tures. TewordtlijO'li£l and Newman and TcwordtEj were 
the first who proposed the complete form of such exten- 
sion. We will use the jaiialogous functional given in a 
more recent publicationcd: 



J^ = a|*| 



2 ' ' 



Kiu^r + 



52 



BHn 



+ 77 1*1 + 



LB'^m^ 



i?'|^|2|D^J-|2 



477 

i?"(V|*|')2 



P D-"* 



Qrot^B. (49) 

The last term rote's was written incil in the equivalent 
form -i rof B(WD* - *D*). 

To account for all the perturbations of the order of t 
one should assume that coefficients u, R' , R" , P, L, and 
Q are temperature independent whereas coefficients a, g 
and K are expanded in t as: 



a = (ai +a2<)t, 

9 = 90+ 9it, 
K = Ko + Kit. 



(50) 



The microscopic BCS values of these coefficients are given 
in Appendix. 

The combination g^l"^ jK that enter in (O) as defined 
by (|6|) parameter k becomes now temperature dependent. 
We keep a notation k = (1 + 7)/V2 for the temperature 
independent part of (ra) and take atj \f2 with 



o' = ffi/50 - 'i.Ki/K^ 



(51) 



as a linear in t contribution to (|6|). Therefore the third 
term in ( p7| ) contains both the perturbation in 7 and in 
t and is written as: 



(7 + ai)(|^|'-l)2 



(52) 



Other perturbation terms of (49) can be substantially 
simplified if one takes into account that they are oper- 
ating with solutions of BJR equation. The final form of 
these terms in dimensionless variables and correspond- 
ing dimensionless coefficients are given in Table I. We 
present also the numerical values of these coefficients cal- 
culated from microscopic BCS theory given in Appendix. 
We comment now, how the perturbation terms were ob- 
tained. 

is rewritten in dimensionless 



1) The term P|D2*| 
units as: 



-i/x \2[{d - ia){d - ia) + [d - ia)(d - ia)]i^\ (53) 
= — i/i |4(9 — ia){d — ia)^^ + 2i[da — 9a)'0| ■ 

Because of (O) and (O) the first term in brackets van- 
ishes and the second term is equal to: 



tiJib^ ivr 



-i^i^i (1 



,|2^2 



y 



(54) 



2) The term _R' |^| |D^| is rewritten in dimensionless 
units as: 

-2tp' IV'I^ (2 |(<9 - m)V'|^ + 2 |(a - m)V'|^) (55) 
= -2tp'{At IV'I^ \(d - m)i/'f - iad \i}\^ + iad \%j:t). 

The first term in brackets vanishes and two other ones 
can be integrated by parts. This leads to: 



tp'\i}}\^ ■2i{da~da)^-tp'b\ 



(56) 



-tp' 1^1^(1 -1^1^). 



3) The term Qrot^B in dimensionless units is rewritten 



as: 



2x2 



trrot^i'Lb) = -iT(V|V'|j 



(57) 



Multiplying (El]) by |V'| and integrating by parts 
find that (V j-^l )^ can be substituted as: 



(v|vi'f-IVl'-|^r, 



(58) 



the same substitution was also done for the term 

i?"(V|^|')2. 

Collecting all the above contributions together and 
omitting the nonessential constant contribution 1 + 7 + 
at — ho we came to the resulting perturbation functional: 

/' = {ho - h,2) IV^I' + (7 - C4t) 1^1' - cet 1^1' , 

(59) 

where parameters Ci are given in Table II. The instability 
field: 



1 + 27 + C2t. 



(60) 



corresponds to the upper critical field that we discuss be- 
low. 

Perturbation functional (p9[) is the principal result of 
the present work. It allows to calculate the properties 
of superconductor with low 7 and select the most sta- 
ble vortex configuration at given 7, ho and t. Although 
functional (59) resembles the extended form of the GL 
functional, it is defined for the restricted set of infinitely 
degenerate vortex solutions |-0(r, ri...r„)| of BJR equa- 
tion (pi|). The important advantage of the functional 
( p9| ) is that, it depends only on few parameters: hc2, C4 
and cg that are the combinations of the coefhcients of 
the extended GL functional (^) as it is given in Table 
II. Moreover, it is not necessary to know the coefficients 
in the starting functional (^) at all. These parameters 
can be considere d as the phenomcnological ones. As will 
be shown in Sec. VII they can be found from experiment. 



The occurring vortex state depends on sign and rela- 
tive strength of the coefficients 7 — C4i and — cg that can 
be both positive and negative since t < and parame- 
ter 7 changes the sign when k goes though l/v^. The 
realistic values of C4 and cq will be discussed in Sec. \Vl% 

(i) When both 7 — c^t and — cgi are positive the mag- 
netic behavior of the superconductor corresponds to the 
generic scenario for a superconductor of type II. The 
dense vortex lattice (^5|) appears continuously from the 
normal state |V'| = at upper critical field hc2 when the 
quadratic term {ho — hc2) IV'I in ( |59| ) becomes unstable. 
The amplitude A^ =< |V'(r)| > of the vortex state and 
the intervortex distance increase with decreasing of the 
applied field. Bel ow th e low critical field hci that will be 
calculated in Sec. VI B all the vortices continuously leave 



the superconductor and the uniform Meissner state with 
\ip\ — 1 becomes stable. 

(ii) When both 7 — c^t and —c^t are negative, the 
functional (p3) corresponds to a type-I superconductor. 
There are only two competing local minima of (|59|): the 
normal state with |^| = and energy 



and the uniform superconducting Meissner state with 
jf/'l = 1 and energy: 



fs = ho - /lc2 + 7 - (C4 + C6)t. 



(62) 



The discontinuous transition between them occurs at the 
thermodynamic critical field he 



he = 1 + -f + {C2 + C4 + CQ)t. 

that is found by equating ( |6l| ) and (|6^). 

(iii) We investigate the case when co effic ients 7 
and 



(63) 



VI] 



c^t 
It will ap- 



-cgi have different signs in Sec 
pear that, depending on the situation, both N-V and V- 
S transition can be either continuous (as in conventional 
superconductors with k ^ 1/ v2) or discontinuous. This 
situation is accessible experimentally either by variation 
of 7 or by variation of t. 



VI. VORTEX STATE: ENERGY AND CRITICAL 
FIELDS 



A. Energy of the vortex lattices and higher critical 

field 



The energy of the regular iV-quanta lattice (p9h can be 
written in terms of the amplitude of the order parameter 

A2 =< \ijf > as: 



/' = {ho - he2)A^ 



-(7-c4i)/34(A)A4-c6i/36(A)A« 



,(JV) 



where /3„ (A) are the structural factors: 



Pk^K^) 



<|^|'>"/2^ 



(64) 



(65) 



that depend both on the amplitude A and on the lattice 
geometry. 

Minimization of (^^ over A gives the complete infor- 
mation about thermodynamic and magnetic properties of 
vortex lattice in a superconductor with k ^ l/v2, pro- 
vided the dependencies Pn (A) are known. We found 
(3n (A) in the whole region of A (0 < A < 1) using 
the numerical solutions of (p^) for square and triangu- 
lar vortex lattices with TV = 1,2 outlined in Sec. IV F. 
The results are shown in Fig.l as function of magnetic 
induction 



6 = 1 - A^ 



(66) 



,(1,2) 



fn=0, 



(61) 



The values of /?„ at A ^ (i.e. in vicinity of hc2) 
are given in Table III. Parameter p\ ' (0) corresponds to 
the parameter (3 introduced in the original publication 
of Abrikosovfl. Close to hd (where A — > 1) the factors 



Pn are tending to 1. The c orresp onding asymptotic 



expression will be given in Sec. VI E| . 

Functional (64) close to hc2 can be interpreted as a 



Landau expansion of the vortex state energy over the 
amplitude A. When 7 > C4,t the quartic in A term is 
positive and the conventional second order transition oc- 
curs at hc2- When 7 < c^t the transition occurs in a 
discontinuous way either to the finite-amplitude vortex 
state or directly to the Meissner state. The concrete re- 
alization of this transition depends on the relative values 
of cq and C4 and will be discussed in Sec. VII . 
Condition: 



7 — C4t = 



(67) 



defines the tricritical point T2 where the discontinuity of 
N-V transition appears. The discontinuity field h*2 is 
large than the critical field hc2- 

Functional (p3) can be alternatively written in terms 
of b as: 

/' = (27 + C2t) iho-l) + {ho - h,2) {h,2 - h) (68) 

+ (7 - Cit)l3i{b){K2-W - C6tMb)ihc2-W- 

Minimization of /' over b gives the induction fo(/io) and 
the lattice energy f'{ha). Comparing the energies of 
square and triangular lattices with TV = 1,2 at given 
^0 (^ci < ho < hc2) we established that one-quanta tri- 
angular vortex lattice always possess the lowest energy. 
However close to hci the energies of different lattices co- 
incide within the calculation accuracy and this conclusion 
becomes less certain. 



B. Diluted lattice, low critical field and vortex 
interaction 

1. Energy of the diluted vortex lattice 

We consider now the diluted lattice of slightly over- 
lapping A^-quanta vortices assuming that the distance I 
between them is much larger then the coherence length 
(i.e. / 3> 1). Such limit usually occurs close to the low 
critical field hd ■ In this approximation the energy of the 
system is written as: 



/' = rs + 



b m 

£n + 



2ttN 



2 2ttN 



UinS) - hob, 



(69) 



where the background energy of the uniform Meissner 
state /g is given by ([64), En is the one- vortex energy, 
b/2'KN is the density of vortices and —hji is the interac- 
tion of the vortex with an external field. The term Umtil) 
represents the interaction between the nearest neighbor 
vortices. The factor m gives the lattice coordination 
number: m = 6 for the triangular lattice and to = 4 for 
the square lattice. The inter-vortex distance / is uniquely 



related with the vortex concentration 2-!tN /b and geom- 
etry of the lattice as: 



lu = (27r7V/6)i/2, l^ = {^TiN/bViY''^. 



(70) 



The type of the V-S transition depends on sign of the 
long-range vortex interaction Uint{l) that, as will be 
shown below, can be both repulsive and attractive. 

When Uint{l) > the situation is the same as for a 
superconductor of type-II: the V-S transition occurs in a 
continuous way at the low critical field hd that is cal- 
culated from the vortex energy sn- The latter can be 
written on the basis of ( p9|) as: 



EN = 27riV[(l + 27 + C2t) - Ci >h - Cit) + Q 'cet]. 



The structural factors for the iV-quanta vortex: 



c 



(N) 



1 

Njo 



(1 - g^j^)rdr 



(71) 



(72) 



are found from integration of the numerical solution 
.gAr(r) of Eq. ( pSJ ) and are given in Table III. Note also 
that: 



AN) 

'52 



2ttN 



= 1. 



(73) 



.(N) 



Factors Cn for large N will be calculated in Sec. VIC 



Lattice factors f3n (p5[) can be expressed via Ca as: 



Pir^ 



1 - d^^b 

(1 - 6)"/2 



(6-.0). 



(74) 



The A^-quanta vortices penetrate into a sample when the 
positive energy en, required for the vortex creation is 
compensated by the negative magnetic contribution —hob 
i.e. above the critical field: 



, (jv) _ £n 

"■cl " 



27riV 



l + (2-Cr)7 



(75) 



+(c2+crc4+crc6)i. 



The low critical field is defined as the lowest field when 
penetration of vortices becomes favorable: 



hci = mm{h[^ '} 



N- 



(76) 



It appears that only the 1-quanta vortices can appear 
in a continuous way since condition of formation of 2- 
quanta vortices written as h^ < hd or as: 



7 > (c4 + 1.89c6)i 



(77) 



is weaker then the condition of continuity of V-S transi- 
tion: Uintil) > 0, derived below (inequahty (|79|)). 

Vortices penetrate inside a superconductor until the re- 
pulsive interaction counterbalances the energy win. The 



penetrated flux is determined by minimization of 
over b that alternatively can be written as: 



f = fs + [/ici - ho 



4ttN 



U,nMb))]b, (78) 



where dependence l{b) is given by (70). 

When Uint{l) < 0, the transition from the Meissner 
phase occurs either to the finite-density vortex state or 
directly to the normal-metal state in a discontinuous way 
that is manifested by the jump of magnetization. The 
detailed scenario of the transition depends on the energy 
balance bet ween these three phases and will be discussed 
in Sec. 



vn 



The situation is simplified however near the tricritical 
point Ti where the long-range part of Uint{l) changes its 
sign from positive to negative. As will be shown in Sec. 
VI B 2 the short-range vortex interaction in this region is 
still repulsive. The minimum of UintiJ) lies at Z 3> 1 and 
one can apply the nearest-neighbor approximation (78). 
The discontinuity field h*i is smaller than hd- 



The given below calculations of the long-range vortex 
interaction are compatible with calculations of the vor- 
tex interaction given intZI in a different way. Consider 
for simplicity the case of two 1-quanta vortices. When 
the distance between vortices is large (/ ^ 1) it is more 
suitable to describe the vortices in terms of slightly over- 
lapping magnetic fluxes produced by these vortices: 



b±^l-gli\v±l/2\ 



(81) 



as was discussed in Sec. IV C 



The vortex energy is provided by the terms \^p\ , \ip\ 



|6 ■ 



and I -01 in the functional ( p9D that can be evaluated as: 

\^P\' = l-b+-b^, (82) 

H\' = {l-b+-b^f (83) 

= (l-6+)2-t-(l-5+)2-l + 26+6_, 

and 



2. Vortex interaction 

The vortex interaction is known to be repulsive in a 
type-II superconductor (7 3> 0) and attractive in a type- 
I superconductor (7 <^ 0). In this Section we calcu- 
late the interaction energy Uint{l) of two A^-quanta vor- 
tices located at ri^2 — ±1/2 at the intermediate values 
of 7. The numerical part of this problem is based on 
solution of the BJR equation ( |2l| ) with the right-hand 
term 27riV5(r — ri) -|-27riV(5(r — r2) and on substitution of 
this solution into the perturbation functional (py). We 
give the analytical treatment of this problem in cases of 
slightly overlapped (Z ^ 1) and of strongly overlapped 
{I ^ 1) vortices. As the result we obtain that vortices 
begin to attract each other at large distances when: 



7 < (c4 -I- 3ce)t. 



(79) 



Below this instability the vortex interaction has the long- 
range attractive and short-range repulsive character and 
vortices form a bounded state. Inequality (|79| ) presents 
also the condition of discontinuity of VS transition. 

With decrease of 7 the equilibrium distance Iq varies 
from infinity to zero. Below another instability point at 



7 < (C4 -H 1.5 



AN) 
S.8 



AN) 
S6 



AN) 
S.6 



AN) 
S4 



C6)t 



(80) 



the interaction is purely attractive and vortices are stuck 
together with formation of 2A-quanta vortex. The re- 
sults about the short-range vortex interaction can not 
be directly applied to study the vortex lattice since the 
nearest-neighbor approximation ( |7q ) is not applicable at 
low vortex separation I . 



iv^r = (1 - 6+ - 6_)3 



(84) 



= (1 - &H 



)=^ + (1 - b+f - 1 + 6b+b- - 3bH+ - 36^&- 



Only 1-01 and I ■01 terms contain the interaction parts: 
&+6_, b'^b- and 6+6^. With the help of ( p6| ) the overlap- 
ping contribution < b^b- > can be estimated with expo- 
nential accuracy as u(Z)e"'" where u{l) is a slow function 
of I. This term is more important at / ^ 1 then decaying 
as e~^' terms < b+b^ >, < b\b^ >. 

Substituting the contributions < 6+6_ > from (83), 
(84) to ( |59| ) one gets the long-range interaction energy 
per one vortex: 



Uintil) = [7 - (c4 + 3c6)t] • u{l)e' 



-41 



(85) 



This interaction is attractive when condition ( [79|) is sat- 
isfied. 

Consider now the short-range part of the vortex in- 
teraction. When vortices are located close to each other 
(Z ^ 1), their order parameter is given by (p4|): 

|^(r)| = 52(r) -t- iff2(r)(lV)2 Ing^ir). (86) 

To calculate the vortex energy one should estimate: 



(87) 



|0(r)|"d2r«27r / g^{r)rdr 



Since g2(j') is an axisymmetric function, the operator 
(IV)^ can be substituted by yV^. Taking finally into 
account the BJR equation V^ In 172 = ffl ^ 1 o^*^ gets: 



(88) 



Tm^. (2) _ (2). 



The interaction energy is calculated on the base of (|5£ 
with respect to the state where vortex cores coincide. 
With help of (|8|) one gets the short-range interaction 
energy: 

U„uil) = [(7 - C4t)(Cf ' - C?') - 1.5c6t(Cf ' - Cf ^)] • T^l\ 

(89) 

This interaction is attractive when condition ( pO| ) is sat- 
isfied. 



VII. H-T PHASE DIAGRAM 

We are now in a stage to discuss the properties of H—T 
diagram of a superconductor with k ~ l/\/2. (We use 
again the dimensiaaal variables.) Partially this question 
was considered irO on the base of Neumann- Tewordt 
extension of GL equations to the low temperature. The 
advantage of our approach is that, it allows to get the 
structure oi H~T diagram in a unified way from a simple 
perturbation functional (p9|). This functional depends on 
three driving parameters: ho, t, and 7 that are controlled 
by experimental conditions and on three phenomcnolog- 
ical parameters: hc2, C4 and cg that can be found from 
experimentm, basing on the relations: 



C. Energy of the normal-superconducting interface 

The profile of the NS interface is given by (^8|) . Bas- 
ing on our perturbation approach we calculate the NS 
interface energy at ho = he as: 

Cns = -(7-C4i)a4 + C6to6, (90) 

where the structural factors a„ are defined as: 



an 



(|^|'-H")dx. 



(91) 



With help of (Ha) these factors can be present in the form 
of definite integrals: 



(e 



lip 



^nxp 



)dx 



V=0 g2<p 



^nip 



1 -e^v 



-dip'^ 



(92) 



zdT], 



/o v^e-2') + 2ry - 1 
and calculated numerically: 

ai ~ 0.55, ae ~ 0.85, ag ~ 1.06. (93) 

One can now express the vortex structural factors (^n 
( [72| ) for iV ^ 1 via the NS interface factors a„. The 
iV-quanta vortex (E3) can be viewed as the cylindric 
normal-state domain of radius rjv ~ \/2N surrounded 
by the NS interface. The domain energy is written as: 



En = TTr^hc + 27rrArcr„s 



(94) 



where irr'j^hc is the energy of the condensate break in- 
side the domain and 2Trr^ans is the domain wall energy. 
Comparison of (p3) and (|7l|) gives: 



c 



(JV) 



1 



N 



(iV>l). 



(95) 



Interesting to note that formula (pq) can be extrapolated 
to small A'^ with an accuracy 5 — 8%. 



— ^ = 1 + 7 - (c4 -(- ce)t, 

-tin 



['dH)^„ 



87r/34(7-C4i)' 



(96) 



(97) 



He 
Hci 



1 + (C4 - 1)7 + [(C4 - 1)C4 + (Ce - l)c6]t, 



(98) 



extracted from (|g), (B3|), (|68|) and (^). The first value 
was also called as Ki(T)7thc second-one as in P 4^(2 K^iT) — 
1) and the third one as 2K3/lnK3l3. 

We extracted the parameters C4 and cg from mag- 
netic measurements in TaNcB and got: C4 ~ 0.30 and 
cg — —0.15. These parameters can be also estimated 
theoretically from the microscopic BCS-eseression for 
coefficients of the extended GL functionalEil. Calcula- 
tions presented in Appendix and in Tables I and II give 
C4 — 0.5, cg = —0.09 for clean superconductor and 
C4 — 0.1, cg = 0.01 for the dirty superconductor. Al- 
though these estimations do not take into account the 
anisotropy of TaN and the electron-phonon retardation 
effects in ECS theory, they give the correct idea about 
the magnitude of coefficients C4 and cg. We assume fur- 
ther that C4 varies from 0.1 to 0.5 and cg from —0.2 to 
0.01; the negative value of cg being more probable. 

The H ~ T diagram of low-7 superconductor at given 
C4 and Cg can be obtained from comparison of relative 
location of critical fields (^0|), (|63|), ( |75[) and character- 



istic critical points that were calculated in SecVl and 



that are resumed in Table IV. To avoid the narrowness 
of the mixed state region and to clearly demonstrate the 
details of the H — T diagram we trace it in the specially 
normalized coordinates T /T^ and H / He where He lin- 
early depends on T: 



He = H'eiTe - T) , 



(99) 



10 



with H'^ = (47r/5)i/2ai/Tc (Q). Topology of the ff - T 
diagram depends on the relative strength of coefficients 
C6 and C4. Three possible scenario of N-S transition can 
be distinguished. 

i) Fig. 2 corresponds to ce,C4 > and to 7 < 0. 
Between Tc and triple point Li defined by condition 
Hc{T) ^ HciiT) or: 



7 = (c4 + 1.75ce)t, 



(100) 



the superconductor behaves like superconductor of type- 
I i.e. the discontinuous N-S transition occurs at critical 
field He- On the left of Li the N-S transition occurs via 
an intermediate vortex state. The V-S transition occurs 
in a continuous way at H = Hd like in type-II supercon- 
ductor. The N-V transition has a discontinuous character 
close to Li. The discontinuity line H*2 terminates in the 
defined by ( pT\ j tricritical point T2 where the forth-order 
term in functional (p4) becomes positive. On the left 
of T2 the N-V transition occurs in a continuous way at 
H — Hc2 and the superconductor behaves like conven- 
tional type-II superconductor. When 7 decreases, points 
Li and T2 are shifted to low temperatures and the super- 
conductor becomes superconductor of type I in a whole 
temperature region. When 7 increases, points Li and 
T2 are shifted to T^. At positive 7 the superconductor 
has a type II behavior. 

ii) Fig. 3 corresponds to C4 > 0, > ce > — C4/3 and to 
7 < 0. Similar to the case i), the vortex state appears on 
the left of the triple point L2 that is defned by condition 
H,2{T)=H,{T),or. 



7 = (c4 + ce)t. 



(101) 



The N-V transition has a continous character and oc- 
cures at H = Hc2- The type of the V-S tranition is 
provided by the location of tricritical point Ti where the 
long-range vortex interaction chages the sign (condition 
ty^^)- The V-S transition is discontinuous between L2 
and Ti at H = H*i and continuous on the left of Ti at 
H = ifci-When 7 decreases the superconductor trans- 
forms to superconductor of type-I whereas when 7 in- 
creases above zero it becomes a superconductor of type- 
II. 

iii) Fig. 4 corresponds to 04 > 0, — C4/3 > cgp-This case 
corresponds to experimental situation in TaNo.When 7 
is slightly less then zero (Fig. 4a), the H — T diagram is 
obtained from the diagram of case ii) by shift of the tri- 
critical point Ti to the region of negative temperatures. 
Like in previous cases, the superconductor transforms to 
a superconductor of type-I with decreasing of 7. When 
7 increases, point L2 goes to Tc and disappears at 7 = 0. 
When 7 becomes positive (Fig. 4b) the V-S transition is 
continuous between Tc and tricritical point Ti that ap- 
pears at Tc at 7 = and is moving to the region of low 
temperatures with increasing of 7. On the left of Ti the 
V-S transition has a discontinous character. At large 7 
the phase diagram transforms to conventional i? — T di- 
agram of superconductor of type-II. 
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APPENDIX: MICROSCOPIC PARAMETERS 

The microscopic BCS parametria of the extended GL 
functional (|4^ ) were calculated intiJ. The first three co- 
efficients do not depend on the purity of superconductor: 



1 ^ 
Tc 



(Al) 



. 7C(3) 1 
16 {-kTY 



i^ 31C(5) 1 
2 64 {-kTY 



(Here v = mpF/^n'^h^ is the density of states). 

Other coefficients depend on the quality of material 
and can be calculated in two limit cases: 



2t2 



K = iy 



R' 



Q = - 



P 



7C(3) v^ti 
48 (7rr)2' 
j^31C(5) v'^h^ 
3 64 (7rr)4' 
J_31C(5)_fiV_ 
35^64C(3) (7rT)2' 



a. Clean limit 

V 31C(5) h\'^ 



i?" 



20 64 (7rT)4 ' 

V 31C(5) v'^h^ 
~ "12 64 (7rT)4' 
__j/31C(5) e^ ft" 



(A2) 



5 64 c2ft2 (ttT)" 



6. Dirty limit 



i. = .^^, P = -.^(44^, (A3) 



48si7rT' 



12-48s2(7rT)2' 



i?' = -V 



Q = - 



2t2 



12-48si(7rr)3' 

1 ft2l,2 



i?" = -ly 



482si(7rT)3' 



807r S1S2 ' 



L = -v 



80 c^h^ sIs2TtT' 



where parameters si and S2 arc the functions of the scat- 
tering times T, Ti, T2 in the s,p and d channels: 



ft ft,l 1, ft,l 1, 

2Ttr 2 T Tl 2 T T2 



(A4) 
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From ( |Al[ ),( |A2| ) and (A3) we calculate the microscopic 
expressions for the defined by (|5l|) parameter a: 



(Jcl = f , 0-d = 0, 

and for the GL parameter n: 



(A5) 



Hcl 



ttTc 



Kd 



3(7C(3))V^ cs, 

(7<(3)l^)V2 |e|;iw2' '^« ^2(^^)1/2 |g|?,„2 

in clean and in dirty superconductors. 
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TABLE I. The terms of the extended GL functional F, their dimensionless counterparts and coefficients in the perturbation 
functional /. Two last columns give the microscopic BCS values of coefficients for clean and dirty superconductors. 



F 



f 



coefficient 



clean limit 



V = 2aiK u/3g 
p' = aiR'/2gK 
p" = aiR"/gK 
fi = aiP/2K^K'^ 
X — ATvaiL/g 
T = 4TTaiQ/2K'^K 



dirty limit 



IL 1*1° 

Qrot^'B 



-tp"|^l'(l-|V^i') 

-a(ivr-ivr) 

-fr|V|''(l- 1-01^) 



V = -■d/2 ! 
p' = -1? ^ -0.45 
p" = -•i?/2 ^ -0.23 
^ = -9i?/5 « -0.82 
A = -9i?/5 « -0.82 
r = -3i9/5 ^ -0.27 



-0.23 

p' = -7rV21C(3) ~ -0.39 
p" = -7rV42C(3) ^ -0.20 
p = -28C(3)/7r'' ^ -0.35 
A = -36T/57rs2 -^ 
r = -12T/57rs2 -> 



i9 = 31C(5)/49C^(3) ^0.454 



TABLE IL Coefficients of the perturbation functional / 
that are collected from the dimensionless terms of Table I, 
their theoretical BCS-values in the clean and dirty supercon- 
ductors and their experimental estimation in TaN 



Ci 






clean 


dirty 


TaN 


C2 = 
C6 = 


= 2cr + p + A 

= -cr + p' + p" + T - 2p - 

^ V — p — p" + ^, - T 


- X 


0.37 
0.50 
-0.09 


-0.35 

0.10 

0.01 


0.30 
-0.15 



TABLE IIL Structural factors /^i'^' of the Af = 1, 2 quanta 
square and triangular vortex lattices close to Hc2 and struc- 
tural factors ("i of one- and two-quanta vortices. 







A = l 






N = 2 






A 




D 


A 




n 


A'^ 


1.16 




1.18 


1.34 




1.43 


4"' 


1.42 




1.50 


1.95 




2.32 






1.58 






1.45 




'>6 




2.00 






1.75 




'>8 




2.34 






1.99 





TABLE IV. Characteristic points at H-T diagram of superconductor with 7 : 



close to zero. 



4th order coefficient in 

energy expansion < at: 

H4T) = Hc2{T) at: 

Hci{T) = Hc2{T) at: 
Short range vortex inte- 
raction is repulsive at 
ans < at 
Hc(T) = H,^(T) at 



7 < Cit 

7 = (C4 + C(,)t 

7= (C4 + I.22c6)i 

7 > (C4 + 1.26C6)t 

7 < (c4 + 1.55c6)t 
7 = (c4 + 1.75c6)t 



r(2) 



Hl['{T) < Hl^tiT) at: 7 > (c4 + 1.89c6)t 



Long range vortex inte- 
raction is attractive at: 



7 < (C4 -I- 3C(i)t 



Tricritical point T2 where at cg > the Hc2{T) 

transition becomes discontinuous 

Triple point L2 where at ce > the Hc{T), 

Hci{T) and Hc2{T) transition lines meet 

Auxiliary point 

Point where 2-quanta vortex decay onto 

two close-lying 1-quanta vortices 

NS interface energy becomes negative 

Triple point L\ where at ce > the Hc{T), 

HciiT) and H*2(T) transition lines meet 

Two separate one-quanta vortices are more 

stable then one two-quanta vortex 

Tricritical point T\ where at ce < the 

Hci{T) transition becomes discontinuous 
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1.2n 
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n ,N=2 




P. 





2.4 




2.2 


n , N=2 /^ 


2.0 


/^ 


1.8 


/^ n , N=1 


1.6 


/^ ^_^^^^^:^^-:----: 


1.4 


//I^-*==^^''^^^^^ A , N=1 . 


1.2 
-1 n 



0.0 0.2 



0.4 _ 0.6 

b 



1.0 



FIG. 1. Stmctm-al factors /Sf ' =< Jt/;!* > / < It/.]^ >^ and 
/3g =< \ip\ >/<|?/'| > for triangle and square N-quanta 
vortex lattices as function of the average magnetic induction 
b. 



X 



1.04- 
1.02- 
1 nn 


Y=-0.05,c,=0.1,Ce=0.01 
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..■■■'''' '-I 
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1 1 1 1 
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T/Tc 



0.8 



1.0 



FIG. 2. H-T phase diagram of superconductor with 
C4,C6 > and with 7 slightly less then zero It includes nor- 
mal (N), vortex (V) and Meissner superconducting (S) phases. 
Solid and dashed lines corresponds to the discontinuous and 
continuous transitions, dotted lines present the auxiliary crit- 
ical fields. Magnetic field is measured in units of the temper- 
ature dependent critical field He ~ H'^iTc — T). 
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FIG. 3. The same as Fig. 2 but for C4 > 0, > ce > -C4/3 
when 7 is slightly less then zero. 
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FIG. 4. The same as Fig. 2 and Fig. 3 but for C4 > 0, 
— C4/3 > Cg when (a) 7 is slightly less then zero (a) and when 
7 is slightly larger then zero (b). 
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